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GAMMA RAY ASTRONOMY BALLOON RESULTS 


C. E. Fichtel, D. A. Kriiffen, and II. B. Ogelman 
NASA/Goddard Space Flight Center, Greenbelt, Md. 


ABSTRACT 

The results obtained from recent high altitude balloon flights of the Goddard 
gamma ray digitized spark chamber telescope are reported and combined with 
earlier work. The data on the galactic center region is reanalyzed in terms of 
a possible line source in the galactic plane. The analysis leads to a value of 
(2.3 ±1.2) x 10" 4 7/(cm 2 sec rad) above 100 MeV for a galactic latitude interval of 
-3° to +3° in the galactic center region; this result in itself does not justify the 
claim of a detected flux, but is consistent with the line intensity of (4.1 ±0.7) x 
10 " 4 yf (cm 2 sec rad) quoted from recent OSO-III results (Clark et al., 1969). The 
flux and energy spectrum of the atmospheric background were measured as a 
function of angle with respect to the vertical. The measured downward flux is in 
reasonable agreement with measurements of other experimenters at several 
energies; however, the upward flux for energies above 100 MeV is a factor of 
about three lower than that measured in the OSO-III experiment mentioned above. 
This difference seems larger than would be expected , but small compared to the 
difference between the galactic gamma-ray flux observed by OSO-III and the 
predicted one. The upper limits to possible point sources of gamma rays are 
summarized and shown generally to fall below the fluxes expected from a 
straight extrapolation of the x-ray spectra. No positive evidence for a point 


source has been found. 


I. INTRODUCTION 


/ 


Current cosmic ray data and radio astronomical measurements indicate that 
relativistic particles constitute a significant fraction of the total energy in the 
environment of various celestial objects and in some regions of our galaxy and 
the universe. Gamma rays are intimately related to relativistic electrons, and 
high energy nuclear processes. Electrons may radiate 7 -rays by the synchro- 
tron process as they are deflected by magnetic fields , by bremsstrahlung as 
they pass through matter, and by the inverse Compton effect when they interact 
with photons. In high energy interactions of baryons, gamma rays may be pro- 
duced as the decay products of hyperons and mesons , such as the pions formed 
in the interaction of cosmic rays with interstellar matter and matter -antimatter 
annihilation. 

The study of high energy photons is necessarily a new field of astrophysics 
since the short path length of these quanta in air make it necessary to place the 
observing instruments near the top of the atmosphere. Further, the fluxes are 
now known to be small, and at present there has been no certain measurement 
of high energy y-rays (E > 30 MeV) from any point source. However, Clark, 
Garmire, and Kraushaar (1968) have obtained positive evidence for a celestial 
7 -ray flux which is anisotropic with a higher intensity in the direction of the 
galactic plane and a maximum in the galactic center region. The flux measured 
in the galactic plane is more than an order cf magnitude greater than that ex- 
pected from the interaction of cosmic rays with interstellar matter, magnetic 
fields , and photons . 
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In this paper we shall summarize gamma ray balloon flight results related 

to several topics including: (a) present balloon flight results on the gamma 

radiation from the atmosphere, coming upwards, from just below the horizon, 

and downward at balloon altitudes , (b) give the upper limit obtained on the flux 

from M-87, and summarize the upper limits from point sources obtained thus 

far, and (c) review the galactic center data reported previously (Fichtel, Cline, 

Ehrmann, Kniffen, and Ross, 1968) in light of the new results of Clark et al. 

* 

(1968). The atmospheric gamma radiation is of interest both from the standpoint 
of providing a means of determining its specific origin and for the in-flight cali- 
bration of gamma-ray satellite experiments. As a background for this discussion 
the detector system will be briefly described , and an account of the method of 
data analysis will be given. The latter is particularly important in terms of 
understanding the upper limits obtained for the possible sources of gamma 
radiation. 


II. EXPERIMENTAL APPROACH 

For high energy gamma rays, the dominant interaction process is the con- 
version to a negatron positron pair. The average angle between the pair tends 
to decrease with increasing energy, because it is primarily determined by elec- 
tron scattering. Information on the direction and energy of the gamma ray must 
be derived from this pair since the penetration power of gamma-rays and their 
low flux relative to other possible masquerading events make shielding difficult 
and generally undesirable. 

A good gamma ray telescope must discriminate against charged particles 
with a high efficiency, give enough information to select gamma rays from other 
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neutral events, permit accurate determination of the gamma ray's arrival direc- 
tion, provide a means of measuring the gamma ray's energy, and, at the same 
time, have a reasonably large detection probability. The detector system used 
in these experiments is shown in Figure 1. The large plastic scintillator anti- 
coincidence dome together with the directional Cerenkov counter is employed 
to restrict the analysis to downward moving particles and discriminate against 
charged particles. The spark chamber satisfies the need for a large volume 
high information content detector to permit selection of the gamma rays and 
measure the properties of the negatron-positron pair. The central plastic 
scintillator together with the Cerenkov counter in coincidence and the plastic 
dome in anticoincidence provides the information to determine whether or not 

the spark chamber should be triggered. The bell shaped anticoincidence dome 
S 

is a solid piece of plastic scintillator, 5/8" thick. 

The spark chamber itself, the associated electronics, and the advantages of 
a digitized spark chamber are discussed in detail in another article (Ehrmann, 

Fichtel, Kniffen, and Ross, 1967); so only the principal features will be given 

♦ 

here. As shown in Figure 1, the spark chamber consists of two units separated 

t 

by the central scintillator. Each unit consists of a series of fifteen plates and 
sixteen 6" by 6" wire -grid modules. Each of these modules is a frame-wire 
assembly consisting of two orthogonal sets of 128 wires each. One orthogonal 
wire grid serves as the high voltage plane of the modular unit, the other as the 
ground plane. The bottom module in each chamber differs from the others in 
that the wires run at an angle of 45° with respect to the others, in order to allow 
the removal of the ambiguity related to which x reading is associated with which 
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y reading when two tracks are present. The spark chamber gas is 88-1/2% 
neon, 10% helium, 1% argon, and 1/2% alcohol and is at a pressure of 1-1/2 
atmospheres. 

Since the spark is formed between the two wire grids in each x-y module, 
there is no spark robbing, as there is when the pair -producing plates serve as 
the electrodes. In the double wire grid system, an x and y reading for each 
spark is obtained from each modular unit with high efficiency even for multiple 
track events, Each of the grid wires in the spark chamber threads a magnetic 
core, which receives and contains its datum of information when the Iiigh voltage 
is pulsed to the chamber plates. The cores are normally in the quiescent, or 
’’reset," state between operations. Spark current in a wire "sets” a core, and 
the "readout" act of determining which cores have been set by the events resets 
the cores. The mean uncertainty in the position of the particle trajectory at any 
level is approximately 0.4 mm. 

Each of the 30 pair producing plates consists of .02 radiation length of gold 
plated unto a .002 radiation length base of aluminum, thus providing a high Z 
plate v^hich is thin enough so that the character of the event is not destroyed by 
absorption and Coulomb scattering of the electrons, and the direction of the inci- 
dent gamma ray can be measured well. At the same time, the large number of 
plates retains a reasonably large detection efficiency. Finally, information on 
the energy of the negatron and positron can be obtained from the multiple coulomb 
scattering in the plates. 

For the flights looking at point sources, the gamma ray telescope is placed 
in a gondola at a fixed angle with respect to the vertical and the gondola itself is 
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oriented with respect to the geomagnetic field to an accuracy of about ±1/2°. 

The orientation is achieved with the aid of two pairs of gas jets supported on 
booms at each end of the gondola. A flux gate magnetometer Is mounted into the 
gondola in a manner so that a zero reading will occur when the gamma ray tele- 
scope axis is in the north-south plane. The signal from the magnetometer is 
used both for orientation data, together with two other magnetometers, and for 
jet activation if a pointing error greater than about a half degree is detected. 

For the flight looking at the earth albedo, the detector was first set in the 
gondola pointing at an angle of 112° with respect to the vertical. After one hour 
exposure at ceiling with this position, the detector was lowered to an angle of 
172° with respect to the vertical. With the given opening angle of the detector, 
sets of angles with respect to the zenith of approximately 95° to 130° and 155° to 
180° could be studied. 


IE. DATA REDUCTION 

Energy, Arrival Direction, and Flux Calculations 

In order to describe the energy measurement, the accuracy of the arrival 
direction, and also the flux calculation, particularly at low energy, it is neces- 
sary to review first the combined effects of multiple coulomb scattering in the 
spark chamber plates and energy loss by the electrons. The theory of multiple 
coulomb scattering has been analyzed in detail by Williams (1939), MoliSre (1947, 
1948, 1955), and others (Goudsmit and Saunderson, 1940; Snyder and Scott, 1949; 
Scott, 1952; Bethe, 1953) and applied successfully for many years to the energy 
determination of charged particles in nuclear emulsions and cloud chambers. 
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More recontly the method 1ms boon applied to multipJato spark chumbors by 
Plnkau (I960, 1967) and extended to overlapping cells by Kniffon (1907). 

Attempts to take energy loss into account have been mado by Kniffon (1967) 
using the average energy and by Pinkau (1968) using an average rate of energy 
loss. For the high nuclear charge material comprising the platos of this oxpori- 
mont radiation losses (bremsstrahlung) predominate over ionization losses above 
10 MeV and ionization losses predominate below 10 MeV. Including the effect of 
scattering, the energy loss tends to be relatively independent of energy below 
about 25 MeV, whereas above that energy it is proportional to the electron's 
energy. Whereas the ionization energy loss is a smooth function, the radiation 
loss is not and can vary significantly in any one case. For tills reason and be- 
cause of the statistical uncertainty of t»ie energy measurement, either the ap- 
proach of Kniffen (1967) or Pinkau (1968) is adequate. 

Figure 2 defines the geometry used in the coulomb scattering analysis, The 
pv value of the electron, for the case when energy loss is neglected, is given 
by the expression (Pinkau, 1967) : 

pv = 16. 5 Q n /< | /3 p J > , (la) 

where 

Q n 2 = ^ d ^A 2 + ) + ■j(n + 1) n(n - 1) A 2 d, (lb) 

and 

< I Q I > - < I x . - 2x + x I > , flcl 

1 ^nn 1 1 ni n (l-l) n (l-3) ' ' 
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whore d ^ 0.0004 cm and 0.75 cm for this spark chambor, n is tho number of 
plates in a coll, p is tho particle momentum, and v is its velocity. Tho in 
this expression is tho result only of coulomb deflections. In practlco there is 
also a noiso signal, called N, which is essentially independent of energy and re- 
sults from fluctuations of tho spark from tho track location and tho minimum 
accuracy for determining x . For this experiment, N was 1.0 mm. Tho meas- 
ured signal 0 is given by the following relationship; 

- <l/ij> 2 + <|N|> a . (2) 

<|N|> may be determined by calibration or by eliminating it using tho overlapping 
cell method (Fowler, 1950); the former method was used here since it (a) permits 
using only n « 1 where that is appropriate and (b) gives a more accurate measure 
of <|n|> on the average. 

In practice, calibration data on electrons of known energy were used to obtain 
experimental curves of <|/3 mn |> vs. electron energy. This was done for each 
deck in the top half of the chamber , where , in each case , all readings arising 
from data taken above that deck were eliminated in determining <|/3 [> and the 
result plotted against the most probable energy ol ine electron at that deck. Thus 
we obtain calibration curves for the determination of the energy of any electron 
formed in the conversion of a gamma-ray at any level in the top half of the 
chamber. In this manner an energy estimate of each electron is obtained (Note; 
since v % c, pv % E), and the energy of the gamma ray is obtained by adding 
them with the uncertainty in the gamma ray energy being the square root of the 
sum of the squares of the uncertainties in the individual energies. The uncertainty 
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in the gamma ray energy typically varies from 30% at the lowest energies to a 
factor of 2 at about 150 Me V. 


The arrival direction of the gamma ray is determined from the knowledge 
of each electron’s energy and direction in the following manner. First, on the 
basis of the energy estimate of the electron, a determination is made of the 
number of plates, m , through which the electron must pass before <|/3 sn J > is 
equal to <|n|> . Then, a least squares fit is made to each electron through the m 
plates, beginning at the apex of the pair. This approach provides the maximum 
information on track location (Kniffen, 1967), unless a point weighting method is 
used (Pinkau, 1968). The latter improves the accuracy somewhat, but greatly 
increases the data reduction time. The projecicad angles of the gamma ray in 
the two orthogonal projections are then calculated as follows is given by the 
expression 


4,r = (“1 4,1 + 4 , 2 ) (3) 

where 4> x . is the projected angle of the jth electrons and to is the statistical 
weight of the <£ xj determination, which is analyzed in Appendix I and given by: 

"1.2 = ct 2 2 i < 4 ) 

where 



and a is the r.m.s. error in the x measurement which is 0.5 mm and k 2 1 is the 
number of plates for which <|/3 g J> = <|n|> . £ d 2 W ni is defined in Appendix I. 
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The accuracy of the determination of the arrival direction can be calculated 
from the lateral displacement function for each electron (Pinkau 1966, 1967). 

The final angular errors quoted are those within which 95% of the true gamma 
ray arrival directions will lie. 

Above about 120 MeV, the arrival direction is determined by the bisector of 
the pair and the uncertainty in arrival direction at very high energies (> 1000 MeV) 
is limited mostly by the track location accuracy. In addition, there is also an 
inherent small angle between the direction of the electron and the primary gamma 
ray due to the pair producing interaction. When both electrons are taken into 
account, the error introduced is about 0.3 of the above errors in each case 
(Stearns, 1949). Since the uncertainties add in a random way, this latter effect 
increases the uncertainty by only about 5%. A curve showing the uncertainty in 
arrival direction as a function of the primary gamma ray energy is shown in 
Figure 3. 

Knowing the balloon gondola angles, the balloon location, and the time, the 
conversion from the angles in the spark chamber to celestial coordinates , or 
angles with respect to the vertical in the case of the albedo and atmospheric 
secondary studies can then be accomplished in a straight forward manner 
(Kniffen, 1967). 

In qrder to determine the absolute value of the flux of gar la rays , it is 
necessary to determine the efficiency for detection as a function of energy and 
the angle of the detector with respect to the direction of interest. At high ener- 
gies, the efficiency is simply a product of the pair production probability and 



the fraction of the sensitive area exposed to the source. At the low energies the 
efficiency is affected both by the energy loss of the electrons in the chamber and 
electron scattering which causes some electrons to be scattered into and out of 
the cone of acceptance. This effect is considered in detail by Kniffen (1967) with 
the net result being to increase slightly the detection efficiency above some 
energy and to cause it to fall to nearly zero at about a gamma ray energy of 
30 MeV for a point source. For a diffuse source such as the atmospheric back- 
ground, the considerations are similar. Figure 4 shows the (area-time-detection 
efficiency) factor as a function of energy for a particular flight and potential 
source. The absolute value depends on the exposure time and the angle of the 
source with respect to the detector axis as a function of time, but the shape of 
the curve is only slightly sensitive to the direction of the source relative to the 
detector axis. 


IV. RESULTS 


(a) Galactic Center Region 

In the introduction it was mentioned that, with an experiment on OSO-in, 
Clark, et al. (1968) detected gamma rays from the galactic plane with a flux 
which was strongest in the region near the galactic center. In an earlier paper, 
results of our group were published which set an upper limit for the flux from a 
point source assumed to be at the galactic center. The balloon flight data in- 
cluding the galactic center region have been reanalyzed in terms of a possible 
line source of finite width centered about the galactic plane. For reference 
Figure 5 shows a distribution of the arrival directions of the observed y-rays 
with measured energies above 100 MeV. Most of these y-rays are, of course, 


10 , 



atmospheric secondaries. Also shown in the figure are the contours of equal 
detecting efficiency and the galactic coordinates. From the figure it is seen that 
a region of the galactic plane from about -10° to -f-25° galactic longitude was 
examined. This is approximately the region of maximum intensity observed by 
Clark et al. (1969), wherein the average line intensity measured was (4.1 ±0.7) 
x 10 “ 4 7 's/(cm 2 sec rad) in this region for energies greater than 100 MeV. 

In the analysis of the balloon flight data obtained in the flight on Dec. 10, 
1966, four different spatial intervals were examined; these were the region 
within the contour for 50% of maximum detection efficiency shown in Figure 5 
from +15° to -15° galactic latitude, ±10° to -10° galactic latitude, +5° to -5° 
galactic latitude, and ±3° to -3° galactic latitude. For each region, the number 
of observed 7 -rays above 100 MeV and above 150 MeV was determined. For the 
energy interval from 30 to 100 MeV, only the first three regions were studied, 
because the angular resolution of the 7 -ray's arrival direction was not adequate 
to justify examining the smallest interval. From the flux calculated in this man- 
ner was subtracted the background gamma-ray flux, which was estimated from 
that observed in the region within the contour for 50% of maximum detection 
efficiency, but excluding the region within ±15° of the galactic center. The re- 
sulting flux was then converted to a line intensity and the results are shown in 
Table I, where the results of a simple subtraction with a one standard deviation 
error are given in Table la and two standard deviation upper limits are given in 
Table lb. The results in Table I by themselves do not justify the claim of a de- 
tected flux, but the line intensity of (2.3 ±1.2) x lo ~ 4 7 's/(cm 2 sec rad) above 
100 MeV for the +3° to -3° interval is certainly consistent with a positive flux 
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Table la 
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and not in disagreement with a line intensity of (4.1 ±0.7) x 10“ 4 7 f s/(cm 2 sec 
rad) observed in the OSO-III gamma ray experiment. 

(b) Atmospheric Background 

Figure 5 shows the flux of gamma rays above 100 MeV measured in this 
experiment as a function of zenith angle. Notice that the intensity rises from 
that at 9 - 0 (straight downward moving gamma rays) to a maximum near the 
horizon primarily because the average amount of material in which cosmic rays 
can produce gamma rays seen by the detector increases approximately as the 
secant of 9 for small angles. A careful calculation must consider gamma ray 
absorption, source function variation with depth in the atmosphere, angular 
dependence, and other factors. The decrease in intensity after 9 passes 90° is 
partly an energy effect. The gamma rays from tt° decay and other sources will 
be traveling at increasingly large angles with respect to the parent cosmic ray 
particles on the average and will therefore have an increasingly large relativistic 
decrease in their energy when their energy is transformed to the observor's 
reference frame. Hence, a larger and larger fraction of the gamma rays will 
have energies below 100 MeV. Also fewer gamma rays are formed along a line 
looking downward than along one looking at the horizon because a larger portion 
of the line is passing through a region with a lower gamma ray source function. 

The changing shape of the energy spectrum as a function of 6 is illustrated 
in Figure 7. In this figure, the y -ray flux has been corrected for the detection 
efficiency as a function of energy, but not for the effect of detector energy 
resolution. The inclusion of the detector resolution affects the shape in only 
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minor respects. Notice that the spectra become steeper as U Increases, 
particularly at low energies. 

The atmospheric gamma ray flux at 3 g/cm 2 was also integrated over the 
whole solid angle to compare It with the lower energy atmospheric photon spectra 
obtained with crystals which inherently integrate over all solid angles. The re- 
sult is shown in Figure 8. It is seen that the gamma ray flux in the region from 
30 to 200 McV lies above the straight line extrapolation; this result is consistent 
with an enhancement expected in this energy interval as the result of decay of 
7r 0, s produced by collisions of cosmic rays with atmospheric nuclei. 

An interesting compaidson can be made between the balloon results discussed 
here and the OSO-ffi results for the upcoming gamma rays from the earth's 
atmosphere. The flux of y-rays coming directly upwards is essentially the 
same at the balloon altitude of this experiment (~ 3 g/cm 2 ) as it is outside the 
atmosphere , basically because the interaction mean free path of both the charged 
cosmic rays and gamma rays is large compared to 3 g/cm 2 . At large angles 
with respect to the downward direction the flux measured at balloon altitude is 
no longer nearly the same as at a satellite altitude, and, of course, there is a 
downward moving flux of gamma rays at balloon altitudes produced by cosmic 
rays interacting in the atmosphere above the balloon. 

The fluxes of upcoming gamma rays (E > 100 MeV) measured on Explorer 
11 (Kraushaar et al., 1965) and on OSO-III (Clark, 1969) at corresponding 
geomagnetic latitudes are (1.9 ±0.5) X 10“ 3 and (10.5 ±1.0) x 10" 3 gamma 
rays/ (cm 2 sr sec) respectively. The flux measured in this work averaged over 


15 



approximately the same angles is (3.7 ±0.8) x 10" 3 gamma rays/ (cm 2 sr sec). 
Notice that the Explorer 11 results fall below the measurement reported here , 

• but the OSO-ni result is appreciably higher. The differences seem larger than 
would be expected on the basis of estimated errors, but are small compared to 
the difference of more than an order of magnitude between the galactic gamma 
ray flux observed by OSO-HI and the predicted one. 

Tlie fluxes measured here may be compared to the downward fluxes measured 
in other balloon gamma ray experiments at similar altitudes. The results are 
displayed in Table II, and the agreement among these experiments is seen to be 
satisfactory, considering that the quoted errors of some of the experiments are 
only statistical. 

(c) Search for a Discrete Source in Virgo 

On February 27, 1968, the detector system was flown in an oriented gondola 
from Palestine, Texas, in a search for possible discrete sources in the direction 
of the Virgo cluster of galaxies. Particular emphasis was placed on examining 
the region about the radio galaxy M-87 , recently observed to be an emitter of 
hard x-rays (Haymes, et al., 1968). 

The celestial distribution of all gamma-rays observed during this exposure 
was determined under the assumption they were extraterrestrial in origin, The 
number of gamma-rays falling within a square angular bin with boundaries defined 
by the 95 percent confidence limits of angular uncertainty in gamma -ray arrival 
direction was determined for gamma-rays with energies greater than 100 MeV 
and 30-100 MeV. 
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*The detectors in all cases sampled a region of the sky within about 3G r from the vertical. 


The observed intensity of atmospheric secondary gamma rays was 5.7 x 
10“ 3 gamma rays/cm 2 sec stcr for 30 < E y < 100 MoV and 6.5 x 10" 3 gamma 
rays/cm 2 soc ster for E y > 100 MoV. The rather high intensity observed for 
measurements at an atmospheric depth of 3 g/cm 2 results from the fact that a 
malfunction in the orientation system caused the detector axis to tilt to an anglo 
of 45° with respect to the vertical. 


For the intensities quoted above, the expected number of gamma-rays in 
the squares of uncertainty are N B (30-100) « 2.0 and N B (> 100) - 0.6. The 
observed number arc? N 0 (30-100) - 3.0 and N 0 (> 100) - 2.0. The 95% confi- 
dence limits on the number of source counts, S, are then given by the relation- 
ship (See Appendix II) . 



(N b + S) N e ~(N a +s) 
N! 



N“0 


n£ c-n d 
N! 


.05 


(5) 


The 95 percent confidence limits obtained are 3.2 X 10" 4 gamma rays/cm 2 sec 
for 30 < E y < 100 MeV and 2.2 x 10" 4 gamma rays/cm 2 sec for E y > 100 MeV. 
These limits , more than an order of magnitude higher than the capability of the 
detector system, resulted from the failure in the orientation system. This 
failure resulted in poor orientation and eventually a total loss of data after only 
1.5 hours of exposure. 


(d) Summary of Upper Limits on Discrete Sources 

The limits obtained from all flights have been reanalyzed in two energy bins 
and are presented in Table III. For comparison, the best existing limits of pre- 
viously reported data are included. 
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Clark ct al. (1968) suggest that possibly a substantial part of the observed 
galactic gamma ray intensity is due to unrosolvcd discroto sources. Ogelman 
(19G9) has subsequently shown that there is reasonably good agreement between 
the observed dependence on galactic longltudo of tho lino intensity of gamma- 
rays > 100 MoV observed by Clark et al. (19G8) and tho distribution of x-ray 
sources in the galactic plane. Table III includes a column containing the gamma- 
ray fluxes expected from the indicated sources if tho observed hard x-ray spectra 
were extrapolated to gamma-ray energies. For those sources for which the 
spectra is not well determined, the x-ray fluxes In the 1 to 10 A region listed in 
the Friedman (1967) survey have been extrapolated using the E" 2,0 relationship 
suggested by Ogelman (1969). No extrapolation is made for the apparently 
thermal source Sco XR-1 nor for sources not known to be x-ray emitters. 
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Discrete Source Limits 
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APPENDIX I 


CALCULATION OF WEIGHTING FACTOR 
FOR GAMMA-RAY ARRIVAL DIRECTION 


Pinkau (1968) has shown that, for a set of k spark readings, x (il recorded 
at levels z n in the spark chamber, the relative probability of x' and <£', the 
initial coordinate and angle with respect to the spark chamber axis, is given by 


P(y',0 


TT* 1 

(x nl “ x' - i A0* ) 2 

/J, ? V 2 a 2 ) ^ 

(£ d 2 2a2 ) 


(A-l) 


where £ d 2 and W n . are complicated expressions involving the composition of the 
scattering plates and the geometry of the detector and a , the root mean square 
reading error in the spark position. For our detector 



2,89 (n* 

(E „) 2 \ 3 



(A-2) 


where n is the number of plates through which the electron has passed between 
readings. For £ d 2 W n . < < 2a 2 , it is clear that P (x', <£') is maximum for a minimum 

k 

V — ' 

value for the numerator of the argument of the exponent > (x n . - x' - iA<£') 2 . 
This is by definition the least squares fit to the spark chamber data points, 

S ivl ”S and ^mox* 

Pinlcau has further shown that, integrating over x', the distribution function 
may now be expressed as 
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(A-3) 


where 


and 


PW ■ *„> « ^ 


* - w ~ O 


a = ]/Qh ^ " * + 2a2) 


where we have neglected energy loss in the plates. 


The distribution for the two electrons of a pair is then given by 


(P r ) (P 2 ) a ex P 




*1 


erf cr n 


where 


S 12 = A(0' 


and this distribution is maximized with respect to $ max for 


d In (p t p 2 ) 

30 * 


= 0 . 


Hence 


(A-4) 


(A-5) 


(A-6) 


(A-7) 


(A-8) 
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f 


yielding 


Mnnx 


_ jL, °± - 1, O'; 


(A-9) 
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APPENDIX H 

CALCULATION OF UPPER LIMITS 


In any balloon exposure in search of a discrete source ofj gamma rays, a 
celestial distribution is obtained from the determination of the arrival direction 
of each observed gamma-ray. Most of the observed quanta result from second- 
ary gamma rays produced by the interaction of energetic charged particle cosmic 
rays with the matter in the atmosphere above the detector. In an angular bin 
surrounding the source, determined by the accuracy with which the gamma ray 
arrival direction can be determined, N 0 counts will be observed. From the 
overall distribution of arrival directions, assumed to be dominated by second- 
aries, and average value N B of atmospheric produced gamma rays would be 
expected in the bin. The probability that a given N 0 and, assumed average source 
strength s will give the number of observed counts N 0 is given by 


P(s) 



e -(N B + «) (N + S )N 0 


N * 


(A- 10) 


The 95% probability that an average source strength s u or less with a distribution 
P (s) will give the observed number of counts is given by: 

No 


fCD 

J»U 


P(s) ds 


£ 


(s, + N b )‘ 


r 


P(s) ds 


No 

Tp «-*■ N> 

j*0 


■= .05 


(A-ll) 


This equation is the same as that used by Hearn (1968). 
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We feel that this equation gives the best possible estimate of the upper 
limit on the average source strength with the available data since there is in- 
sufficient information to determine whether fluctuations from the expected 
values are due to variations in the source contribution or the background 


contribution. 
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FIGURE CAPTIONS 

Figure 1. Schomatic drawing of digitized spark chamber gamma-ray telescope. 

Figure 2. Scattering coordinates and parameters in the spark chamber. 

Figure 3. Uncertainty in gamma ray arrival direction - 95% confidence as a 
function of energy. 

Figure 4. Area-time-detection efficiency factor as a function of energy. 

Figure 5. Distribution of the arrival directions of the observed gamma-rays 
with measured energies above 100 MeV. The dark lines represent 
the indicated galactic latitude, and the galactic longitude is indicated 
by marks along the 0° galactic latitude line. The dashed curves are 
contours of equal area-solid angel-collection time and the numbers 
by these curves indicate the percentage of the maximum area-solid 
angle-collection time for this flight. 

Figure 6. Gamma-ray flux at a balloon altitude of 3 g/cm 2 as a function of the 
angle with respect to the vertical. 

Figure 7. Energy spectrum of gamma rays observed at a balloon altitude of 

3 g/cm 2 as a function of the angle with respect to the vertical. 

■ • ’ ‘ ’ *■ - - , * . ....... , 

Figure 8. Energy spectrum of atmospheric background gamma-ray flux inte- 



grated over all solid angles. 
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